Entanglement, quantum phase transition and scaling in XXZ chain 
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Motivated by recent development in quantum entanglement, we study relations among concur- 
rence C, SU<j(2) algebra, quantum phase transition and correlation length at the zero temperature 
for the XXZ chain. We find that at the SU(2) point, the ground state possess the maximum concur- 
rence. When the anisotropic parameter A is deformed, however, its value decreases. Its dependence 
on A scales as C — Co — C\(A — l) 2 in the XY metallic phase and near the critical point (i.e. 
1 < A < 1.3) of the Ising-like insulating phase. We also study the dependence of C on the cor- 
relation length £, and show that it satisfies C — Co — l/2£ near the critical point. For different 
size of the system, we show that there exists a universal scaling function of C with respect to the 
correlation length £. 
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Quantum entanglement, as one of the most intrigu- 
ing feature of quantum theory, has been a subject of 
much study in recent years, mostly because its nonlocal 
connotation^ is regarded as a valuable resource in quan- 
tum communication and information processing 0, 0- 
For example, an entangled state, such as a singlet state 
(I Tl) — I IT))) can be used for the realization of 
teleportationQ. On the other hand, as with other re- 
sources, such as free energy and information, one would 
like to know how it can be quantified and controlled. 
For the first problem, much efforts have been devoted 
to develop a quantitative theory of entanglement, in- 
clude entanglement of formation^, ||| 0, @ , which is re- 
garded as its basic measure. For the second problem, 
many authors H E3, EH E3, E3, E3, EM EH E3l tried to 
built bridge between quantum entanglement and physi- 
cal models by inve stig ating their entanglement in both 
the ground state |oL Eg) and thermal state fl^. fl5| . 

Very recently, the intriguing issue of the relation be- 
tween entanglement and quantum phase transition have 
been addressed [TH ITsj . For a spin- 1/2 ferromagnetic 
chain, Ostcrloh et. al. , reported that the entanglement 
shows scaling behavior in the vicinity of quantum phase 
transition point EH as induced by a transverse magnetic 
field. Vidal et. al. tried to establish a connection be- 
tween quantum information and condensed matter the- 
ory by studying the behavior of critical entanglement in 
spin systems. So it is believed that the entanglement 
of the ground state, like the conductivity in the Mott- 
insulator transition 20| and quantum Hall effect, and 
magnetization in the external-field-induced phase tran- 
sition, is also plays a crucial role to the understanding 
of quantum phase transition. On the other hand, group 
theory as well as symmetry of the system are parts of 
the foundation of quantum mechanics |U|22|, the knowl- 
edge of its presence often makes it easy to understand the 



physics. Thus the study of entanglement at the ground 
state and its relation to the group theory will not only 
have a contribution to experimental realization, but also 
enrichs our physical intuition of quantum theory. 

The main focus of the present paper is to study the 
properties of ground state concurrence of an antiferro- 
magnetic XXZ chain. We show that the competition be- 
tween quantum fluctuation and ordering will lead to max- 
imum value of concurrence at the isotropic point. This 
observation could also be clarified from the point view 
of g-deformation theory. The concurrence's dependence 
on anisotropic parameter A is presented both numerical- 
lly and analytically. The relation of the concurrence to 
the correlation length £ in the Ising-like insulating phase, 
as well as the scaling behavior around the critical point 
A = 1 where the Metal-insulator quantum phase transi- 
tion occurs, are also discussed. Thus our result not only 
manifest interesting physical phenomenon, but also es- 
tablish non-trivial connection between the quantities in 
quantum information theory and critical phenomenon, 
correlation length in condensed matter physics and quan- 
tum group theory [22. 

The Hamiltonian of the XXZ chain with periodic 
boundary conditions reads 
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H(A) = £[ ff f af +1 + a?cr? +1 + Aafaf +1 ] 
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where N is the number of sites, a a (a = x, y, z) are Pauli 
matrices, and A is a dimensionless parameter character- 
izing anisotropic interaction. The Hamiltonian is invari- 
ant under translation, therefore, the entanglement be- 
tween arbitrary two neighbor sites is a uniform function 
of site index. At A = 1, Eq. (JIJ has SU(2) symmetry. 
While A 7^ 1, it becomes g-deformed SU(2) algebra with 
A = (q+q~ 1 )/2. Together with the Z 2 symmetry, we can 
have [H, S z ] = 0, which result in that the reduced density 
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matrix pm+i) of two neighbor sites is of the form[j| 



Pi(i+i) 



fu+ 
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(2) 



in the standard basis | TT)> I t-L), I IT) J ID- Since the en- 
ergy of a single pair in the system is E/N — tr[pm + ijHi], 
where Hi is the part of Hamiltonian between site I and 
I + 1, due to the translational invariance. Considering the 
definition of entanglement, we can easily obtain that the 
concurrence of XXZ chain can be calculated as^J 113 



C=- max[0, \E/N - AGf +1) | - Gfg +1) - 1]. (3) 

where G zz l+1 ^ is the correlation function. So we not only 
need to know the energy of the system, but also the be- 
havior of correlation function. 

It is well known that the present model can be ex- 
actly solved by Quantum Inverse Method 0j H3 , and 
its energy spectra are determined by a set of spin rapidi- 
ties Ai, A2, . . • , Am, which describe the kinetic behavior 
of a state with M down spins. They are the solution of 
Bethe-ansatz equation 



/ sinh7(Aj + i) 
\s'mhj(Xj — i) 
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M 
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sinh7(Aj — A; + 2i) 
sinh7(Aj — A; — 2i) 



(4) 



where the parameter 7 arises from the anisotropic scale 
A, i.e., A = cos 27. The regime < A < 1 is character- 
ized by real positive 7 while the regime 1 < A by pure 
imaginary 7 with positive imaginary part. When 7 — > 0, 
the above secular equations reduce to the well known one 
for isotropic Heisenberg model. 

Taking the logarithm of the above equation, we can 
have a set of transcendental equations for {Xj}, in which 
the energy level is determined by a set of quantum num- 
ber {Ij}. For the ground state, {Ij} are consecutive 
integer or half-odd-integer centering around zero, and 
M = N/2. Then the ground state energy of the sys- 
tem can be calculated either by solving the Bethe ansatz 
equations numerically for finite size system, or by solving 
integral equation of density function of A in the thermo- 
dynamic limit. Once the A dependent eigenenergy E(A) 
is obtained, the correlation function is simply the first 
derivative of E(A)/N with respect to A. 

For the XXZ model, there exist two different phases 
at the ground state, i.e., metallic phase: < A < I 
and insulating phase: A > 1, which is resulted from 
that the former is gapless while the later is gapful. The 
critical point of quantum phase transition locates at the 
isotropic point A = 1 at which the concurrence is just 
a simple function of ground state energy per sites, i.e. 
0.386. If we regard a z as a 'coordinate', then the first 
two terms in Eq. represents the 'kinetic' energy caus- 
ing the quantum fluctuations of a z , and the last term 



represents the 'potential energy' that causes the order- 
ing of a z . In the Ising limit A — ► 00, the ground state 
has the Neel long-range-order, which results in that the 
concurrence is zero. When A becomes smaller but still 
large than 1, the quantum fluctuation plays more and 
more important role, then the Neel state is no longer an 
eigenstate of the Hamiltonian. This fluctuation between 
two neighboring sites enhances the value of off-diagonal 
term z in their reduced density matrix pm+\), then the 
entanglement becomes larger and larger. On the other 
hand, at the free particle(XX) limit where A = 0, the 
spin-flip term dominates the system completely, and all 
spins flip freely on lattice sites. For a certain site j, the 
probability of spin up and down is the same, regardless 
the spin state of its neighbor. Thus the state | ||) will 
not lower the energy, but share the same probability with 
I ||) or I ||). This phenomenon will result in a rela- 
tively large u + or u~ in the reduced density matrix of 
two neighbor sites, as well as a relative smaller C. On 
the contrary, once the anisotropic interaction 1 > A > 
is turned on, the value of u + and u~ is lowered. So the 
concurrence is enhanced. Hence the competition of quan- 
tum fluctuation and ordering must results in a maximum 
concurrence at a certain point. Comparing with the ori- 
gin of metal-insulator transition in the present model, 
which also arise from the competition of fluctuation and 
ordering, it is natural to infer that the point we want 
here is just the isotropic point, i.e. A = 1, as illustrated 
in the Fig. ^ This case is very similar to the forma- 
tion of Kondo effect, in which the competition between 
spin singlet formation and thermal conductivity leads to 
a minimum conductivity at the Kondo temperature. The 
idea can also be applied to the entanglement of arbitrary 
two sites, such as the concurrence C; m between site I 
and m. Only when the competition between their in- 
teraction and fluctuation reaches a counterbalance, the 
concurrence C\ m reaches its maximum. 

From the quantum group theory point of view, at A = 
1 point, the ground state is SU(2) singlet in which the two 
neighboring sites try to form antisymmetric pair, as (| || 
) — I IT))/v2. In the g-deformed region, the Hamiltonian 
(|T|) can be rewritten in terms of Temperly-Lieb operators 



H = NA 



N 

2£T, 



(5) 



where Tjj+i = {—q \ 1; 1, -q} in the basis | ||), | ||). 
Define q-deformed antisymmetric state \<p q ) = (| || 
) — l\ IT))/\/l + l 7 2 i then the operator Tj j+i can be ex- 
pressed as Tjj + i = — — \4>q) {4>q I- If A > 1, the lowest 
energy state favor the formation of g-deformed antisym- 
metric state between two neighboring sites |2fj, unlike the 
case of A = 1 where it favors antisymmetric state, which 
obviously leads to the decrease of concurrence between 
two neighboring sites. When the deformation parameter 
q becomes very large, it tends to the Neel state. On the 



other hand, the 



breaks the local translational invari- 



ance, from the point view of spinless fermions model, the 



formation of \<j> q ) develops charge-density-state (CDW) 
at the ground state, which is gapped and low symmetric. 

We show the concurrence as a function of A in Fig. 
n which is obtained by solving both the Bethe ansatz 
equations for 1280 sites system numerically, and the in- 
tegral equation for infinite length system (We obtained 
the same result). As we expect, the ground state at 
the isotropic point possesses the maximum concurrence. 
Thus symmetry of the Hamiltonian plays a central role 
in determining the concurrence of its ground state. And 
the trend of curve can be easily understood based on the 
above argument. On the other hand, a challenge and 
non-trial problem is to quantify the concurrence around 
the critical point. In the XY metallic phase and near 
the critical point (i.e. 1 < A < 1.3) of the Ising-likc 
insulating phase, it is amazing that C can be described 
by 

C = Co-d(A-l) 2 , (6) 

very well, where 

C = 21n2-l ~ 0.386, 

d = 21n2-i---4- 0.047, (7) 

as illustrated in the inset of Fig. ^ Hence around A = 1, 
the critical exponents of the anisotropic term is 2. As we 
know, the present model can be transformed into spinlcss 
fermions model by Jordan- Wigncr transformation. For 
the free particle case, it is easy to obtain that the ground 
state energy and the correlation function G zz are 4/-7T 
and 4/-7T 2 respectively. In the large A limit, we find the 
concurrence scales like C tx 1/A. One can also express 
C in terms of deformation factor q via the relation q = 
A ± VA 2 - 1. It has the form 

C = C -^{q^-q-^f. (8) 

around the critical point. In XY metallic phase, it we 
define q = e^, it becomes 

C-C , -4C 1 sin 4 ^. (9) 

Now we study the scaling behavior in the Ising-like 
insulating phase by considering the correlation length. 
Though the scaling study of metal-insulator transition 
based on the analysis of spin stiffness has proposed 
recently [53], and though everyone believe there must ex- 
ist some relation between correlation and entanglement, 
the scaling of concurrence, and its dependence on the 
correlation length still remains an open and interesting 
problem. By analyzing the finite chain system, one can 
obtain the correlation length as a function of A in an 
easy wav |27ll28| . It has the form 

1 /e = 7 + ^^^tanh(2n 7 ), (10) 

n—l 
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FIG. 1: Representation of concurrence C as a function of 
A, obtained by solving three sets of Bethe ansatz equa- 
tions of N — 1280 sites system numerically. It is clear that 
the concurrence reach is maximum at the critical point A = 1. 




FIG. 2: Representation of concurrence C as a function of 
Here £ is in unit of lattice constant. 

In the A — > 1 limit, it has a good approximation, as l/£ oc 
(A — l) 2 . Clearly, the correlation length is independent 
of system size, its behavior is shown in the inset of Fig. 
121 The dependence of C on £ is represented in Fig. [3 in 
which the solid line is obtained by solving Bethe ansatz 
equations for 1280 sites system numerically. For the value 
of £ bigger than 4, i.e. l/£ < 0.25, there exist a simple 
relation between C and £, which scales 

C = C -^. (11) 

The above equation imply that the concurrence of does 
not have a long-range effect, in another way, we can say 
that a smallish system, such as TV = 20, can well describe 
the behavior of concurrence of large system, as illustrated 
in Fig. |3| Compared with the scaling of spin stiffness |27|. 
the present one is more perfect, that is the concurrence is 
almost independent of the system size when L > 10. So 
we can conclude that for finite size system, there exist a 
scaling function, which is independent of L and scales like 
Eq. (|llfl in large £ limit. Only when L < 10, the finite 
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FIG. 3: Representation of concurrence C as a function of £ 
for different system sizes. 

size effect becomes very clear (See the inset of Fig. OJl. 
Moreover, for small system, concurrence in even number 
sites and odd one is different. The former is usually larger 
than the later due to the frustration effect happens in 
odd sites system with periodic boundary condition. For 
example, for 3 sites system, the two singlet formations 
between sites 1, 2 and between sites 2, 3 breaks singlet 
formation of sites 3, 1. When L becomes large, this effects 
can be neglected and the concurrence in two case is the 
same. 



In summary, we have investigated the ground state 
concurrence of the XXZ chain. We pointed out that 
the competition between quantum fluctuation and Neel 
ordering will lead to a maximum value of concurrence 
at the isotropic Heisenberg point. Based on the Bethe 
ansatz solution, we exactly obtained the dependence of 
C on the parameter A in a wide range around the critical 
point, and obtained numerical result in the whole range. 
We established the relation between the concurrence and 
deformation factor q of quantum group in the Ising-like 
insulating phase. It is now clear that (/-deformed per- 
mutation generator favors the formation of a deformed 
ground state, which has a relatively smaller concurrence. 
Moreover, the relation between the concurrence and the 
correlation length was studied both numerically and an- 
alytically. We found that there exists a universal scaling 
behavior for finite(not small) size system, and it satis- 
fies a simple relation C oc l/2£ in the region close to the 
critic point. 
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